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Abstract 

Collinear factorized perturbative quantum chrome-dynamics (pQCD) model predictions are 
compared for p + Pb at 5. 02 A TeV to test nuclear shadowing of parton distribution at the 
Large Hadron Collider (LHC). The pseudorapidity distribution the nuclear modification factor 
(NMF), RpPb(y = 0,Pt < 20 GeV/c) = drippy/ {N co \\(b)dn pp ) and the pseudorapidity asymmetry 
Y asym.(PT) = Rp P b(PT,V < 0) / ' Rp Pb (PT , V > °) are computed using HIJING/BB v2.0 model and 
a pQCD improved parton model kTpQCD_v2.0 which embedded generalized parton distribution 
functions (PDFs). These results are updated calculations of those presented in Ref. [l|. 

PACS numbers: 12.38.Mh, 24.85. +p, 25.75.-q, 24.10Lx 
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I. INTRODUCTION 



In this note we show predictions for moderate px < 20 GeV/c observables in p + Pb 
collisions at 5.02 A TeV at the Large Hadron Collider (LHC) usingthe HIJING/BB v2.0 



model 



3| and a pQCD improved parton model kTpQCD_v2.0 [4], [5J which embedded 
generalized parton distribution functions (PDFs). All model details are extensively discussed 
in the literature and we focus only on the updated results for pseudorapidity distributions, 
nuclear modification factor, R p pb(T),PT,b) = dnpPb/(N co a(b)dn P p) and the pseudorapidity 
asymmetry, Y^ sym {p T ) = Rp Pb (p T ,r] < 0) / Rp Pb (p T , r) > 0). These predictions are testable 
with a short 5.02A TeV run (10 6 events). 

II. NUCLEAR SHADOWING AND JET QUENCHING AT LHC ENERGIES 

Monte Carlo models as HIJING1.0 fl, HIJING2.0 Q and HIJING/BB2.0 1-3] have been 
developed to study hadron productions in p + p, p + A and A + A collisions. They are es- 
sentially two-component models, which describe the production of hard parton jets and the 
soft interaction between nucleon remnants. The hard jets production is calculated employ- 
ing collinear factorized multiple minijet within pQCD. A cut-off scale po in the transverse 
momentum of the final jet production has to be introduced below which (px < po) the 
interaction is considered nonperturbative and is characterized by a finite soft parton cross 
section a so f t . Jet cross sections depend on the parton distribution functions (PDFs) that are 
parametrized from a global fit to data M, M. 

nn 

Nucleons remnants interact via soft gluon exchanges described by the string models [91.1101] 
and constrained from lower energy e + e, e + p, p + p data. The produced hard jet pairs and 
the two excited remnants are treated as independent strings, which fragments to resonances 
that decay to final hadrons. Longitudinal beam jet string fragmentations strongly depend 
on the values used for string tensions that control quark-anti-quark (qq) and diquark-anti- 
diquark (qqqq) pair creation rates and strangeness suppression factors (7 S ). In the HIJING1.0 
and HIJING2.0 models a constant (vacuum value) for the effective value of string tension is 
used, kq = 1.0 GeV/fm. At high initial energy density the novel nuclear physics is due to 
the possibility of multiple longitudinal flux tube overlapping leading to strong longitudinal 
color field (SCF) effects. Strong Color Field (SCF) effects are modeled in HIJING/BB2.0 
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by varying the effective string tensions value. SCF also modify the fragmentation processes 
resulting in an increase of (strange)baryons which play an important role in the description of 
the baryon/meson anomaly. In order to describe p+p and central Pb + Pb collisions data at 
the LHC we have shown that an energy and mass dependence of the mean value of the string 
tension should be taken into account Moreover, to better describe the baryon/meson 
anomaly seen in data a specific implementation of JJ loops, has to be introduced. For 
a detailed discussion see Refs. j^, [J. Similar results can be obtained by including extra 
diquark-antidiquark pair production channels from strong coherent fields formed in heavy- 



ion collisions 



All HIJING type models implement nuclear effects such as nuclear modification of the 
partons distribution functions, i.e., shadowing and jet quenching via a medium induced 
parton splitting process (collisional energy loss is neglected) pj. In the HIJING1.0 and 



HIJING/BB2.0 models Duke-Owen (DO) parametrization of PDFs [12j is used to calculate 
the jet production cross section with pt > Po- In both models using a constant cut-off 
Po = 2 GeV/ c and a soft parton cross section cr so ft = 54 mb fit the experimental p + p data. 
However, for A + A collisions in HIJING/BB2.0 model we introduced an energy and mass 

nn 

dependence of the cut-off parameter, p (s,A) [2|, |3] at RHIC and at the LHC energies, in 
order not to violate the geometrical limit for the total number of minijets per unit transverse 
area. 



In HIJING2.0 Ml model that is also a modified version of HIJING1.0 



6] the Gluck-Reya- 



a 

Vogt (GRV) parametrization of PDFs [3] is implemented. The gluon distributions in this 
different parametrization are much higher than the DO parametrization at small x. In 
addition, an energy-dependent cut-off po(s) and a so f t (s) are also assumed in order to better 
describe the Pb + Pb collisions data at the LHC. 

One of the main uncertainty in calculating charged particle multiplicity density in Pb + 
Pb collisions is the nuclear modification of parton distribution functions, especially gluon 
distributions at small x. In HIJING type models one assume that the parton distributions in 
a nucleus (with atomic number A and charge number Z), f a /A{x,Q 2 ), are factorizable into 
parton distributions in a nucleon ( f a /N) an d the parton(a) shadowing factor (S a /A), 

f a /A(x, Q 2 ) = S a/A (x, Q 2 )Af a/N (x, Q 2 ) (1) 

In our calculations we will assume that the shadowing effect for gluons and quarks is the 
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same, and neglect also the QCD evolution (Q 2 of the shadowing effect). At this stage, the 
experimental data unfortunately can not fully determine the A dependence of the shadowing. 
We will follow the A dependence as proposed in Ref . {(J and use the following parametrization, 

c / \ fa/A(x) 

= 1 + 1.191og 1/6 A[x 3 - 1.2x 2 + 0.21x] 

-^ 1/3 -m-^^V^e-^\ (2) 

Sa = 0.1, (3) 

The term proportional to s a in Eq. |2]determines the shadowing for x < Xq, (where Xo=0.1) 
with the most important nuclear dependence, while the rest gives the overall nuclear effect 
on the structure function in x > xq with some very slow A dependence. This parametrization 
can fit the overall nuclear effect on the quark structure function in the small and medium 



x region j(5]. Because the remaining of Eq. [2] has a very slow A dependence, we will only 
consider the impact parameter dependence of s a . After all, most of the jet productions occur 
in the small x region where shadowing is important: 

s a {b) = s a h -{l-b 2 /R\\ (4) 



where Ra is the radius of the nucleus, and s a = s q = s g = 0.1 The LHC data {3] indicate 
that such quark(gluon) shadowing is required to fit the centrality dependence of the central 
charged particle multiplicity density in Pb + Pb collisions. This constrain on quark(gluon) 
shadowing is indirect and model dependent. Therefore, it is important to study directly 
quark(gluon) shadowing in p + A collisons at the LHC. In contrast, in HIJING2.0 a 
different A parametrization ((A 1 / 3 — l) - 6 ) and much stronger impact parameter dependence 
of the gluon (s g = 0.22 — 0.23) and quark (s q = 0.1) shadowing factor is used in order to 
fit the LHC data. Due to this stronger gluon shadowing the jet quenching effect has to be 
neglected 

Note, all HIJING type models assume a scale- independent form of shadowing parametriza- 
tion (fixed Q 2 ). This approximation could breakdown at very large scale due to dominance of 
gluon emission dictated by the DGLAP [l^| evolution equation. At Q = 2.0 and 4.3 GeV/ c, 
which are typical scales for mini-jet production at RHIC and LHC respectively, it was shown 



that the gluon shadowing varies by approximately 13% in EPS09 parametrizations 15] 
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FIG. 1: (Color online) (a) HIJING/BB2.0 predictions of charged particles pseudorapidity distribu- 
tion (dN c h/ 'drj) for minimum bias (MB) p+Pb collisions at 5.02^4 TeV. Solid curves includes fixed 
Q 2 shadowing functions from HIJING1.0 6j and SCF effects, while the dashed curve has SCF 
effects but no shadowing. The dotdashed curve are the results withouth SCF and Shadowing, (b) 
Ratio RpPb(rj) calculated assuming iV co n(MB) = 6.4 

III. HIJING/BB PREDICTIONS 

Figure [1] shows HIJING/BB2.0 predictions of the global observables dN^/drj and -R p pb(f?) 
= (dNpp h ) I 'drj) I (N co )idNpp / drj) characteristics of minimum bias p + Pb collisions at 5.02 A 
TeV. The predictions for p + p are also shown. Minijet cutoff and string tension parameters 
Po = 3.1 GeV/c and k = 2.0 GeV/fm for n + Pb are determined from fits to p + p and A + A 

u n 

systematics from RHIC to the LHC (see Refs. [2j, yj, for details). Note, these calculations 
assume no jet quenching. 

The absolute normalization of dN^/dr] is however sensitive to the low px S, 2 GeVj_c 



nonperturbative hadronization dynamics that is performed via LUND 9| string JETSET 10] 
fragmentation as constrained from lower energy e + e, e+p, p+p data. The default HIJING1.0 
parametrization of the fixed Ql = 2 GeV 2 shadow function leads to substantial reduction 
(solid histograms) of the global multiplicity at the LHC. It is important to emphasize that 
the no shadowing results (dashed curves) are substantially reduced in HIJING/BB2.0 relative 
to no shadowing prediction with default HIJING/1.0 from Ref. 6|], because both the default 
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minijet cut-off po = 2 GeV/c and the default vacuum string tension kq = 1 GeV/fm (used 
in HIJING1.0) are generalized to vary monotonically with centre of mass (cm) energy per 
nucleon y/s and atomic number, A. As discussed in 2|, [3] systematics of p + p and Pb+Pb 
multiparticle production from RHIC to the LHC are used to fix the energy (y/s) and the A 
dependence to a cut-off parameter po(s, A) = 0.416 y/s*' 191 A 0A28 GeV/ c and a mean value 
of the string tension k(s,A) = Ko (s/sO) 04 A 0167 GeV/fm j^j]. The above formulae lead to 
Po = 3.1 GeV/ c and k = 2.1 GeV/fm at 5.02A TeV for p + Pb collisions. For p + p collisions 
at 5.02 TeV we use a constant cut-off parameter po pp = 2 GeV/c and a string tension value 
of K, pp = 1.9 GeV/fm. 
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FIG. 2: (Color online) (a) Minimum bias transverse momentum distributions at mid-pseudorapidity 
\r]\ < 0.8 predicted by HIJING/BB2.0 with (solid histogram) and without (dashed histogram) 



HIJING1.0 shadowing functions 



a. 



The results for p + p collisions at 5.02 TeV (dotted histogram) 
are also included, (b) The mid-pseudorapidity nuclear modification factor of charged hadrons R p pb 
from HIJING/BB2.0 model. The solid and dashed histograms have the same meaning as in part 
(a). 

Note, even in the case of no shadowing shown in Fig. [TJ the increase to po — 3.1 GeV/c 
from p = 2 GeV/ c (value used in p + p at 5.02 TeV) causes a significant reduction by a 
factor of roughly two of the minijet cross section and hence final pion multiplicity. This 
reduction of minijet production is also required to fit the low charged particle multiplicity 
growth in A + A collisions from RHIC to LHC (a factor of 2.2) 0. 
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We interpret this as additional phenomenological evidence for gluon saturation physics 
not encoded in leading twist shadow functions. The px > 5 GeV/c minijets tails are un- 
affected but the bulk low < 5 GeV/c multiplicity distribution is sensitive to this extra 
energy (\/s) and A dependence of the minijet shower suppression effect. It is difficult to 
relate po to saturation scale Q sa t directly, because in HIJING hadronization proceeds through 
longitudinal field string fragmentation. The energy (a/s) and A dependence of the string 
tension value arises from strong color field (color rope) effects not considered in CGC phe- 
nomenology that assumes kx factorized gluon fusion hadronization. HIJING hadronization of 



minijets is not via independent fragmentation functions as in PYTHIA [10j, but via string 
fragmentation with gluon minijets represented as kinks in the strings. The interplay be- 
tween longitudinal string fragmentation dynamics and minijets is a nonperturbative feature 
of HIJING type models. The approximate triangular (or trapezoidal) rapidity asymmetry 
seen in the ratio R p pb{v) sloping downwards from the nuclear beam fragmentation region at 
negative pseudorapidity 77 < —5 toward l/N co a in the proton fragmentation region (77 > 5) 
is a basic Glauber geometric effect first explained in Refs. 13, ll8[ and realized via string 
fragmentation in HIJING. 



T3 



80 
70 
60 
50 
40 
30 
20 
10 
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FIG. 3: (Color online) (a) HIJING/BB2.0 predictions of charged particles pseudorapidity distribu- 
tion (dNch/drj) for central 0-20 % p+Pb collisions at 5.02A TeV. The solid, dashed and the dotted 
histograms have the same meaning as in Fig. [TJ 

In Fig. [2] are displayed the predicted transverse spectra and nuclear modification factor 
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for charged hadrons at mid-pseudorapidity, \rj\ < 0.8. Including shadowing and SCF effects 
reduces R p j>b from unity to about 0.7 in the interestin g 5 to 10 GeV/c region close to the 
prediction of Color Glass Condensate model (KKT04) 19(. A similar nuclear modification 



factor is found 20j using leading order (LO) pQCD collinear factorization with HIJING2.0 
parameterization of shadowing functions [8|, GRV parton distribution functions (nPDF) 
from Ref. 13|, and hadron fragmentation functions from Ref. 21]. 

Figure[3]and Fig.HJshow HIJING/BB2.0 predictions of the global observables dN^/drj and 
R P Pb(v) = (dNpp h ) I drf) I (N co \\dNpp / drf) characteristics of central (0-20 %) p + Pb collisions 
at 5.02v4 TeV. The number of binary collisions iV co ii ~ 12. Figure [3] and Fig. H]The results 
are similar with those presented in Fig. [T] and Fig. [2J In a scenario with SCF and shadowing 
effects the RpPbiPr) are slightly lower that in MB event selection. In contrast in a scenario 
with SCF efects and without shadowing the model predict no suppression for high pt > 5 
GeV/c charged particles ( i.e., ~ 1). 
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FIG. 4: (Color online) (a) Minimum bias transverse momentum distributions at mid-pseudorapidity 
\r]\ < 0.8 (Fig. H^) and nuclear modification factor for central 0-20 % p + Pb collisions at 5.02 ^4TeV 
predicted by HIJING/BB2.0 model. The solid and dashed histograms have the same meaning as 
in Fig. [2J 
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IV. CALCULATIONS WITH kT P QCD_v2.0 MODEL 



ogically enhanced, perturbative QCD 
% |5j. The main feature of this model 



The kTpQCD_v2.0 code is based on a phenomeno 
improved parton model described in details in Refs. 
is the phenomenologically generalized parton distribution function in order to deal with 
the non-perturbative effects at relatively low-x (i.e. small p?) values. The model includes 
the so called intrinsic-/cT parameter as phenomenological corrections for non-perturbative 
effects determined by data from wide energy range of nucleon-nucleon (mainly pp) collisions. 
Moreover, within the framework of this model, the broadening of the intrinsic-/^ in proton- 
nucleus (pA) or nucleus-nucleus (AA) collisions is related to the nuclear multiple scattering. 
This can gen erate enhancement of the nuclear modification factor, the so called Cronin 



effect 



22 



23[, which appears within the 3 GeV/c < < 9 GeV/c region from SPS to 



RHIC energies. 

The kTpQCD_v2.0 code calculates the invariant cross section for hadron production in 
pp, pA or AA collisions, which can be described at LO or NLO levels in the /^-enhanced 



pQCD-improved parton model on the basis of the 
Monte Carlo based integration of the convolution 

rl-q-W'c dv r 
<VW/z c V(l-V) Jvw/vz c W 



771 h 

^T3 — 

d 6 PT 



abc 



? actorization theorem. The code provides 



x 



Q 2 ) fb/ P (xb, k T 6, Q* 
— d[l-w) H K ab , c (s,v,w,Q,Q R ,Q F ) 

dV 7T 



D 



c K^ci 



Q 



(5) 



where we introduced the 3-dimensional generalized parton distribution functions in a fac- 
torized form, 

/(x,k T ,Q 2 ) = f(x,Q 2 )-g(k T ) . (6) 

Here, the function f(x,Q 2 ) represents the standard 1-dimensional LO or NLO PDF as a 
function of momentum fraction of the incoming parton x at factorization scale Q, da/dv rep- 
resents the Born cross section of the partonic subprocess ab — > cd, K abtC (s, v, w, Q, Qr, Qf) 
is the corresponding higher order correction term, and the LO or NLO fragmentation func- 
tion (FF), D^(z c , Q F ), gives the probability for parton c to fragment into a hadron, h with 
momentum fraction z c at fragmentation scale Qf- We use the conventional proton level 
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(S, V, W) and parton level (s, v, w) kinematical variables of next-to-leading order calcula- 



tions ( for more details see Refs. 



24 



2*. 



In this analysis we consider fixed scales: the factorization and the renormalization scales 
are connected to the momentum of the intermediate jet, Q = Qr = n-p q (where p q = Pt/z c ), 
while the fragmentation scale is connected to the final hadron momentum, Qp = n-p?. The 
value of k = 2/3. 

We introduce the phenomenologically generalized 3-dimensional (3-D) PDF which is as- 
sumed to be factorized to a standard PDF (at its Q scale) and a 2-D initial transverse- 
momentum distribution, g(kr) of partons containing its 'intrinsic-/^' parameter as in 
Refs. 0, [5I, Q, Q- We demonstrated the success of such a treatment at LO level in Ref. 
4j, and a Kj et -based NLO calculations in Refs. 28|, |29|. In our phenomenological approach 
the transverse-momentum distribution is described by a Gaussian, 



<7(k 



7l(k 2 T ) 



(7) 



Here, (k^) is the 2-D width of the kp distribution and it is related to the magnitude of the 
average transverse momentum of a parton as (k^) = A(kx) 2 /iT. In order to reproduce the 
nucleon-nucleon collisions at relatively low-x, we assumme (kV) = 2.5 GeV 2 /c 2 . 

As a standard 1-D PDFs, MRST( at central gluon (eg)) 301 was used in eq. (EJ). For 
FF we use the most recent parameterizations KKP from Ref. 31[ . These sets of PDF and 
FF can be applied down to very small scales (Q 2 ~ 1.25 GeV 2 ). Therefore, we can perform 
calculations at relatively small transverse momenta, p T > 2 GeV for our fixed scales. 



A. Initial State Nuclear Effects in pA and AA Collisions 

Proton-nucleus and nucleus-nucleus collisions can be described by including collision ge- 
ometry, saturation in nucleon-nucleon (NN) collision number, and shadowing inside the 
nucleus. In the framework of Glauber picture, the cross section of hadron production in 
nucleus-nucleus collision can be written as an integral over impact parameter, h: 

Eh A <L = / d^ dV tA(r) M |„ _ r |) . E , j<(<4U(4U) (8) 
d p T J dp 

where pp cross section on the right hand side represents the cross section from eq. (EJ), but 
with an increased widths compared to the original transverse-momentum distributions ([7]) 
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in pp collisions, as a consequence of nuclear multiscattering (see eq. (jHJ)). Here tA.(b) = 
J dz pA{b, z) is the nuclear thickness function (in terms of the density distribution of nucleus 
A, pa), normalized as J d 2 btA(b) = A. For the Pb nuclei Woods - Saxon formula was applied. 

The initial state broadening of the incoming parton's distribution function is accounted 
for by an increase in the width of the Gaussian parton transverse momentum distribution 
in eq. (|7j): 

(k 2 T ) pA = {k 2 T } pp + C ■ h pA (b) (9) 

Here, (k 2 -,)^ is the width of the transverse momentum distribution of partons in pp collisions, 
h p A{b) describes the number of effective NN collisions at impact parameter b, which impart 
an average transverse momentum squared C. The effectivity function h p A{b) can be written 
in terms of the number of collisions suffered by the incoming proton in the target nucleus, 
va(P) = cr NN tA(b), where <tnn is the inelastic NN cross section: 



h pA (b) 



u A (b) - 1 u A (b) < v m 
v m — 1 otherwise 



(10) 



We have found that for realistic nuclei the maximum number of semihard collisions is 3 < 
v m < 4 with C = 0.4 GeV 2 /c 2 . 

Furthermore, the PDFs are modified in the nuclear environment by the 'shadowing' 
effect jsl, (32)- 34]. This effect and isospin asymmetry are taken into account on average 
using a scale independent parameterization of the shadowing function S a /A{ x ) adopted from 
Ref. bd: 



S a /A{%) 



Z fa/ P (X,Q 2 ) + ( 1 " ~ A ) fa/n(x,Q 2 " 



(11) 



where f a / n (x, Q 2 ) is the standard 1-dimensional PDF for the neutron and Z is the number o 
protons. In the present work, we display results obtained with the EKS99 [32[, EPS08 [33 1 
and HKN 34 j parameterization, and with the updated HIJING parameterization The 
former three has an anti-shadowing feature, while the latter one incorporates different quark 
and gluon shadowing, and has an impact-parameter dependent and an impact-parameter 
independent version. The impact-parameter dependence is taken into account by a term 
oc (1 — b 2 /R\), which re- weighs the shadowing effect inside the nucleus. 
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FIG. 5: (Color online) Predictions updated at 5.02A TeV from Ref 



. HQ 



for central - 20% 



(b < 3.5 fm) pPb at midrapidity. Predictions with with DGLAP Q 2 evolved HIJING 



EPS08 



33(], and HKN 



from PHENIX 



38 



39] 



EKS99 



341 ] shadowing are also shown. The data for d + Au collisions at 200 GeV are 



B. Results on Nuclear Effects 



The nuclear modification factor, Rpp b (pT) is presented on Fig [5] for charge-averaged 



hadrons, using various shadowing functions, S a / Pb (x). On this plot HIJING 



EPS08 



EKS99 



32|, 



33], and HKN 34J shadowing parameterizations are plotted. For the case of HIJING 



parameterization both impact-parameter dependent /independent versions are plotted. Cal- 
culations were made for yfs = 5.02 ylTeV, |?7| < 0.3, and — 20% central pPb. Note, due to 
the lack of 6-dependence of the above mentioned shadowing functions, this is equal to the 
minimum bias results as well. All details and parameters are the same as in Ref. ij. See 



more details in Refs. 



35 



-w.\. 



We perform calculations on rapidity asymmetry defined as: 



y (pt) 

asym Kf-L J 



E h d z o h p Pb/d p T \ v<0 _ R% Pb (p T ,r)< 0) 



(12) 



E h d 3 <T$Pb/d>,\ v>0 R h pPb (p T , V > 0) ' 
For these calculations we used kTpQCD_v2.0 with various type of shadowing functions. 



12 




kTpQCD v2.0 with 
b-indep. HIJING 2.0 

EKS99 



' b-indep 
Shadowing * \ 
+ Multiscattering ' 



HKN shadowing 

MinBias HIJING BB 2.0 
w/o Shadowing 
w SCF 4 w JJ 



MinBias pPb-> h LHC 5.02 ATeV at rj= 0.3-0.8 

5 10 15 20 25 30 

Pt (GeV/c) 



0.6 



kTpQCD v2.0 with 
0-20% b-dep. HIJING2.0 
0-20% b-indep. HIJING 2.C 
0-20% EKS 



0-20% HKN shadowing 



•'0-20% b-indeffW 
HIJING Shadowing 
+ Multiscattering 



0-20% HIJING BB2.0 
w/o Shadowing 
w SCF + w JJ 



0-20% pPb^ h LHC 5.02 ATeV at tj= 0.3-0.8 

5 10 15 20 25 30 

Pt (GeV/c) 



FIG. 6: (Color online) Predictions for m {px) updated at 5.02 ATeV from Refs 



asym \ 



3. Q,Q 



for 



MB events selection (left panel) and central O — 20% pPb collisions (right panel). Compared are 



fixed Q 2 deeply shadowed HIJING predictions 
Predictions with DGLAP Q 2 evolved EKS99 



81 with and without impact parameter dependence. 



31 EPS08 



331 ] shadowing are shown. The results 



obtained using HIJING/BB v2.0 model (stars) are also included. 



Calculations are at = 5.02ATeV using the 0.3 < \rj\ < 0.8 rapidity range to the back- 



ward and forward directions. Results are an extension of those published in Ref. [37| on 
effects close-to-midrapidity region. The results are plotted on Figure M for unidentified 
charged hadron production in minimum bias (MB) pPb collisions (left panel) . One com par e 



the ca 



culations within kTpQCD_v2.0 code including HIJING EKS99 |32|, EPS08 |33|, 



HKN 34] shadowing parameterizations. The results using HIJING/BB v2.0 mode 



are ob- 

SCF and J J loops effects [sj. Due 
to soft physics embedded in the model the predictions are different compared with pQCD 
inspired model kTpQCD_v2.0. The reason for this is under study now and will be presented 
elsewhere. Similar results are plotted for central (0 — 20%) p+Pb collisions in the right 
panel. 
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V. CONCLUSION 



In conclusion, even with a small sample of 10 6 events the study of RpPbiPr) or central 
relative to peripheral NMF (-Rcp(Pt)) could provide a definitive constraint on nuclear shad- 
owing implemented within different pQCD inspired models and CGC saturation models, 
with high impact on the interpretation or reinterpretation of the bulk and hard probes for 
nucleus-nucleus (Pb+Pb) collisions at LHC energies. 
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